Multiple Superconducting Phases in New Heavy Fermion Superconductor PrOs4St>i2 
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The superconducting gap structure of recently discovered heavy fermion superconductor 
PrOs4Sbi2 was investigated by using thermal transport measurements in magnetic field rotated 
relative to the crystal axes. We demonstrate that a novel change in the symmetry of the supercon- 
ducting gap function occurs deep inside the superconducting state, giving a clear indication of the 
presence of two distinct superconducting phases with twofold and fourfold symmetries. We infer 
that the gap functions in both phases have a point node singularity, in contrast to the familiar line 
node singularity observed in almost all unconventional superconductors. 

PACS numbers: 74.20.Rp, 74.25. Fy, 74.25.Jb, 74.70. Tx 



In almost all superconducting (SC) materials known 
to date, once the energy gap in the spectrum of elec- 
trons opens at the SC transition, only its overall ampli- 
tude, and not the shape and symmetry around the Fermi 
surface, changes in the SC phase §. The vast major- 
ity of superconductors have conventional s-wave pairing 
symmetry with an isotropic SC gap which is indepen- 
dent of directions over the entire Fermi surface. Over 
last two decades, unconventional superconductivity with 
gap symmetries other than the s-wave has been found in 
several classes of materials, such as heavy fermion (HF) 
compounds^, high-T c cuprates||, ruthenateQ, and or- 
ganic compounds^. There, the strong Coulomb re- 
pulsion leads to a notable many-body effect and often 
gives rise to Cooper pair states with angular momen- 
tum greater than zero. Unconventional superconductiv- 
ity is characterized by anisotropic superconducting gap 
functions belonging to nontrivial representation of the 
crystal symmetry group which may have zeros (nodes) 
along certain directions. The nodal structure is closely 
related to the pairing interaction of the electrons. It is 
widely believed that the presence of nodes in almost all 
unconventional superconductors discovered until now are 
a signature of magnetically mediated interaction, instead 
of conventional electron-phonon mediated interaction. A 
common feature in the HF superconductors discovered 
until now is that they all have line nodes in the gap func- 
tions, parallel or perpendicular to the basal planes. 

Recently a new HF superconductor PrOs4Sbi2 with 
filled skutterudite structure has been discovered Q. This 
material should be distinguished from the other uncon- 
ventional superconductors, in that it has a non-magnetic 
ground state of the localized /-electrons in the crystalline 
electric field (most likely doublet T 3 state) || fj], || . The 
origin of HF behavior in this compound (m* ~ 50m e , m e 
is the free electron mass) was discussed in the light of the 
interaction of the electric quadrupole moments of Pr 3+ , 
rather than local magnetic moments as in the other HF 
superconductors, with the conduction electrons. There- 



fore the relation between the superconductivity and the 
orbital fluctuation of /-electron state (i.e. quadrupole 
fluctuation) has aroused great interest; PrOs4Sbi2 is a 
candidate for the first superconductor mediated neither 
by electron-phonon nor magnetic interactions. Hence it 
is of the utmost importance to determined the symme- 
try of the SC gap. Very recently heat capacity jump at 
T c in zero field was reported to be broad even in high 
quality samples 0. This behavior was attributed to the 
double superconducting transition, implying a possible 
novel feature of the superconducting state. 

Here we studied the SC gap structure of PrOs4Sbi2 by 
the thermal transport measurements. Cooper pairs do 
not carry entropy and therefore do not contribute to the 
thermal transport. As a result, thermal conductivity is a 
powerful probe of the low energy excitations of unpaired 
quasiparticles (QPs). Especially the measurements of the 
angular variation of the thermal transport in a magnetic 
field rotated relative to the crystalline axes proved to be 
a uniquely powerful probe for determi ning the location of 
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port on what seems to be an exception to the well-known 
observation concerning the robustness of the symmetry 
of the SC gap function and the presence of line nodes, 
providing compelling evidence for the novelty of the SC 
state of PrOs4Sbi2- 

Our PrOs4Sbi2 single crystals were grown by the Sb- 
flux method and had T c =1.82 K with a very sharp resis- 
tive transition AT C /T C < 0.01. Clear de Haas van Alphen 
oscillations were observed at 0.4 K, confirming the high 
quality of the sample |^2f . In our samples, the heat ca- 
pacity jump at T c in zero field showed a structure of some 
sort |2j| , similar to Ref. j?J . We measured the c-axis ther- 
mal conductivity k zz (the heat current q|| c) on the sam- 
ple with a rectangular shape (0.40x0.37x2.20 mm 3 ) by 
the steady-state method. To apply H with high accu- 
racy relative to the crystal axes, we used a system with 
two superconducting magnets generating H in two mu- 
tually orthogonal directions and a 3 He cryostat set on 
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FIG. 1: Main panel: T-dependence of the c-axis thermal con- 
ductivity k zz in zero field and above H C 2- Insets; Upper: 
T-dependence of k zz (H — Q)/k zz (H = 2.5 T) near T c . Lower 
left: H-dependence of k zz ( £Z"||[100]) at 0.4 K. Lower right: 
schematic figure of the H=H (sin 9 cos <fi, sin 9 sin <f>, cos 9) ro- 
tation. While rotating <f>, 9 was kept constant. 



a mechanical rotating stage at the top of the Dewar. 
By computer-controlling the two magnets and rotating 
stage, we were able to rotate H with a misalignment of 
less than 0.5° from each axis, which we confirmed by the 
simultaneous measurements of the resistivity. 

Figure 1 shows the T-dependence of k zz in zero field 
and above H C 2 (~ 2.2 T at T— K). In this tempera- 
ture region, the electronic contribution to k zz dominates 
the phonon contribution. In the upper inset of Fig. 1, 
Kzz(H — T) normalized by k zz (H = 2.5 T) is plotted. 
Upon entering the SC state, a double shoulder structure 
is seen at T ~1.6 K and 1 K, which may be related to 
the double SC transitions j?| g3| . The lower left inset of 
Fig. 1 shows the if-dependence of k zz . n zz is nearly lin- 
ear in H over almost entire field range below H C 2- This 
behavior is in a marked contrast to the exponentially 
slow increase of the thermal conductivity with field ob- 
served in s-wave superconductors at H <C H C 2 ■ The steep 
increase of k zz in PrOs4Sbi2 is a strong indication that 
the thermal transport is governed by the delocalized QPs 
arising from the gap nodes p4|. 

A strong additional evidence for the existence of the 
nodes is provided by the angular variation of k zz in H 
rotated within the a6-plane. The theoretical understand- 
ing of the heat transport for superconductors with nodes 
has made significant advances during past few years p5[ . 
The most significant effect on the thermal transport for 
such superconductors comes from the Doppler shift of the 
QP energy spectrum (e(p) — > e(p) — v s ■ p) in the circu- 
lating supercurrent flow v s p6| . This effect becomes im- 
portant at such positions where the local energy gap be- 
comes smaller than the Doppler shift term (|A| < |v a -p|), 
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FIG. 2: (a)(b) Angular variation of k zz in H rotating within 
the afc-plane (9 — 90° in the lower right inset of Fig. 1) at 
0.52 K above and below H c2 {~ 2.0 T). 



which can be realized in the case of superconductors with 
nodes. The maximal magnitude of the Doppler shift at a 
particular point strongly depends on the angle between 
the node direction and H . For instance, when H is ro- 
tated within the basal plane in quasi two dimensional 
d-wave superconductors, the Doppler shift gives rise to 
the fourfold oscillation of the density of states (DOS). 
In this case, the DOS shows the maximum (minimum) 
when H is applied to the antinodal (nodal) directions 
Jl6| , [nj. This fourfold pattern is also present in the 
transport properties, and thus the thermal conductivity 
can be a powerful probe of the gap structure. In fact, 
a clear fourfold modulation of the thermal conductivity, 
which reflects the angular position of nodes, has been 
observed in high-T c cuprate YBa 2 Cu 3 7 ||Hl |ll], HF 
CeCoIn 5 @, organic k-(ET) 2 Cu(NCS) 2 |13), and boro- 
carbide YNi 2 B 2 C fllij , while it is absent in the fully 
gapped s-wave superconductors fll| . Here we apply this 
technique to PrOs4Sbi2. 

In the present experiments, k zz was measured by rotat- 
ing H= H(sin9 coscf>, sin sin 0, cos#) conically around 
the c-axis as a function of 0, keeping 9 constant |Q. Here 
9=(q,H) is the polar angle and (f> is the azimuthal angle 
measured from the a-axis (see lower right inset of Fig. 1). 
Figures 2 (a) and (b) display the angular variation of 
k zz (H,4>) in H rotating within the a6-plane (9 = 90°) 
at T=0.52 K. The measurements have been done in ro- 
tating H after field cooling at <f> = —90°. The open 
circles show k zz (H, <j>) at H=1.2 T and 0.5 T which are 
obtained under the field cooling condition at each angle. 
Values of K zz (H,<j>) obtained by different procedures of 
field cooling are nearly identical, indicating that the field 
trapping effect related to the vortex pinning is negligi- 
bly small. Above H c2 (~2.0 T at 0.5 K) K zz (H,<f>) is 
essentially independent of <fi. A clear fourfold variation is 
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FIG. 3: (a)The amplitude of twofold (open circles) and four- 
fold (filled circles) symmetries plotted as a function of H/H C 2 
at T=0.52 K. For details, see the text. (b)Angular variation 
of n zz (H, 4>) at 6 = 90°, 45° and 30° measured by rotating H 
( \[ibH\=1.2 T) conically around the c-axis as a function of 
4> (see the lower right inset of Fig. 1). 



observed just below H c2 down to H ~ 0.8 T. However fur- 
ther decrease of H below 0.8 T causes a rapid decrease of 
the amplitude of the fourfold symmetry and eventually 
a discernible fourfold symmetry is not observed below 
0.7 T. At the same time, the twofold symmetry grows 
rapidly. 

Figure 3(a) shows the ^/-dependence of the ampli- 
tude of twofold and fourfold symmetries. The amplitude 
are obtained by decomposing k zz (H, (j)) into three terms 
with different symmetries; n zz {H,4>) = kq + «20 + k^, 
where kq is a (^-independent term, K24, — C24, cos 2<fi, and 
= C40 cos are terms with twofold and fourfold 
symmetry with respect to ^-rotation. We attempted to 
fit k zz (H,4>) using functions which possess the two and 
fourfold symmetries other than cos 2(f) and cos 4<^ and 
found that the independence of both amplitudes are in- 
sensitive to the choice of the functions. It is apparent that 
the transition from the fourfold to twofold symmetry in 
^-rotation occurs in a narrow field range at H/H C 2 ~ 
0.4, which is deep inside the SC state. Both symmetries 
coexist in a narrow field range, possibly due to the inho- 
mogeneity. 

It should be emphasized that the fourfold anisotropy 
of the Fermi velocity vf, which is inherent to the cu- 
bic band structure of PrOs4Sbi2, is quite unlikely to be 
an origin of the observed twofold and fourfold symme- 
tries because of the following reasons. First, the fact 
that k zz (H, (j)) is ^-independent above H C 2 strongly indi- 
cates that anisotropy of the Fermi surface does not pro- 
duce a discernible anisotropy in the thermal transport. 
Second, the observed twofold symmetry at lower field is 
lower than the symmetry of the Fermi surface. Third, the 
amplitude of the anisotropy of H C 2 within the afo-plane, 
which is less than 1%, is too small to explain the observed 



fourfold variation of n zz (H,(f)) below H C 2- Fourth, the 
calculations based on the Kubo formula indicate that the 
leading term in the thermal conductivity is related to the 
anisotropy of the gap function and anisotropy of vf will 
only enter as a secondary effect. These consideration lead 
us to conclude that the observed symmetries originate 
from the gap nodes. Therefore the change from fourfold 
to twofold symmetry provides a direct evidence for the 
change of the gap symmetry. The fact that K zz (H,(f>) 
possesses minima at <\> = ±90° and 0° at high field phase 
and at <f> = ±90° for low field phase immediately leads 
to the conclusion that the nodes are located along the 
[100]- and [010]-directions in the high field phase, while 
they are located only along the [010]-direction in the low 
field phase. 

We discuss here why the lower twofold symmetry ap- 
pears in the cubic lattice with higher symmetry. The 
twofold symmetry indicates the absence of a double do- 
main structure, suggesting a mechanism which fixes a 
single domain. As shown in Fig. 2(b), the same twofold 
symmetry is observed in K zz (H,(j)) measured under the 
field cooling condition at each angle. This fact excludes 
the possibility in which the magnetic field fixes the do- 
main. Another mechanism, such as one dimensional dis- 
locations produced in the process of the crystal growth, 
may be a possible origin for the single domain structure. 
Further study is required to clarify its origin. 

Having established the presence of nodes, the next 
question is their classification. As demonstrated in 
Ref. jl4|, |TJ , the angular variation of k zz can distinguish 
between the point and line nodes, by rotating H conically 
around the c-axis with a tilted angle from the afo-plane. 
In this case the amplitude of the fourfold and twofold 
symmetries decrease with decreasing 9 for point nodes, 
while they essentially maintain the same magnitude for 
the line nodes. This can be explained by noting that 
rotating H tilted from the afe-plane does not point to 
the nodes for the point node case, while rotating H at 
any tilt angle always cross the nodes for the line node 
case. Figure 3(b) displays the angular variation of k zz 
at fixed 9. The amplitude at 9 = 45° and 30° is smaller 
than that at 9 = 90°. Similar results were obtained for 
the twofold symmetry. These results are in favor of point 
nodes. We note that these results are consistent with the 
recent NQR experiments which exclude the presence of 
the line nodes p7|. 

We next discuss the nodal structure inferred from the 
present results. Unfortunately whether the point nodes 
are present along the [001]-direction remains unresolved. 
This is because the measurements of function 
of the angle 9 rotated across the c-axis is necessary to 
identify the nodes along the [001]-direction. In such mea- 
surements, however, the dominant twofold oscillation in 
k zz (9), arising from the difference in transport lifetime 
of QPs traveling parallel to and normal to the vortices 
p6| , p7| , pj[ p0[ , obscures the part of the signal coming 
from the nodal structure along the [001]-direction. We 
shall therefore be content here with a discussion of the 
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FIG. 4: The phase diagram of the superconducting gap sym- 
metry determined by the present experiments. The filled cir- 
cles represent the magnetic field H* at which the transition 
from fourfold to twofold symmetry takes place. The open cir- 
cles represent H C 2- The area of the gap function with fourfold 
symmetry is shown by pink (^1-phase) and the area of the gap 
function with twofold symmetry is shown by blue (B-phase). 

nodal structure based on a group theoretical considera- 
tion. It is unlikely that the SC gap function has four point 
nodes in the cubic Th crystal symmetry without taking 
into account the accidental mixing between the different 
basis functions in the SC wave function, regardless of spin 
singlet or triplet symmetry. This naively suggests that 
the gap function at high field phase has six point nodes, 
while that at low field phase has two point nodes. More 
detailed study is necessary to clarify the nodal structure 



along [001]-direction. 

The H — T phase diagram of the SC symmetry deter- 
mined by the present experiments is displayed in Fig. 4. 
The filled circles represent the magnetic field H* at which 
the transition from fourfold to twofold symmetry takes 
place. Near T c the amplitudes of both symmetries were 
strongly reduced but were sufficient to resolve if*. The 
il*-line which separates two SC phases (high field A- 
phase and low field -B-phase) lies deep inside the SC state. 
This line seems to terminate at temperature slightly be- 
low T c . It is tempting to consider that the end point co- 
incides with lower transition temperature of the double 
transitions discussed in Ref. 0. Its clarification moti- 
vates further investigations. The only example of a super- 
conductor with multiple phases of different gap symmetry 
so far has been UPt3 [p8| . However, the gap functions 
of UPt3, in which the interaction leading to supercon- 
ductivity is magnetic and the antiferromagntic symme- 
try breaking field reduces the symmetry from hexagonal 
to orthorhombic, are known to be very complicated. On 
the other hand, the SC gap functions of PrOs 4 Sbi 2 are 
expected to be much simpler owing to the cubic lattice 
symmetry. 

In summary, we investigated the gap structure of heavy 
fermion superconductor PrOs4Sbi2 using thermal trans- 
port measurements. Two distinct SC phases with dif- 
ferent symmetries, the phase transition between them, 
and the presence of point nodes-all these highly unusual 
nature highlight the unconventional superconductivity in 
PrOs4Sbi2 and open up a new realm for the study of the 
superconductivity with multiple phases. They also im- 
pose significant restrictions on possible theoretical mod- 
els, including those with a novel mechanism of neither 
phonon nor magnetic origin of pairing. 
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